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Key role of orbital degree of freedom to understand the magnetic structure of uranium compounds 
is discussed from a microscopic viewpoint by focusing on typical examples such as UMGas (M=Ni 
and Pt) and the mother compound UGa3. By analyzing an orbital degenerate Hubbard model 
constructed from the j-j coupling scheme in the cubic system, we obtain the phase diagram includ- 
ing several kinds of magnetic states. Especially, in the parameter region corresponding to actual 
uranium compounds, the phase diagram suggests successive transitions among paramagnetic, mag- 
netic metallic, and insulating Neel states, consistent with the experimental results for AuCu3-type 
uranium compounds. Furthermore, taking account of tetragonal effects such as level splitting and 
reduction of hopping amplitude along the z-axis, an orbital-based scenario is proposed to understand 
the change in the magnetic structure from G- to A-type antiferromagnetic phases, experimentally 
observed in UNiGas and UPtGa 5 . 

PACS numbers: 71.27. +a, 75.30.Kz, 75.50.Ee, 71.10.-w 



I. INTRODUCTION 

Recently, /-electron compounds with the HoCoGa 5 - 
type "115" tetragonal crystal structure have been inten- 
sively investigated in both the experimental and theoret- 
ical research fields of condensed matter physics. Among 
such compounds, interesting magnetic properties have 
been reported for UMGas , where M is a transition metal 



ion 
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In particular, neutron scattering experi- 
ments have revealed that UNiGas exhibits the G-type an- 
tiferromagnetic (AF) phase, while UPdGas and UPtGa 5 
have the A-type AF state^S Note that G-type indicates 
a three-dimensional Neel state, while A-type denotes a 
layered AF structure in which spins align ferromagneti- 
cally (FM) in the ab plane and AF along the c axis^ It is 
quite interesting that the magnetic structure is different 
for U-115 compounds which differ only by the substitu- 
tion of transition metal ions. 

U-115 materials have also been found to differ in their 
magnetic anisotropy. For M=Ni, Pd, and Pt, both \a 
and \c exhibit Curie- Weiss behavior, but Xa is somewhat 
larger than Xcr 1 where Xa and Xc are magnetic suscep- 
tibilities for magnetic field parallel to the a and c axes, 
respectively. This anisotropy increases in the sequence 
UNiGa 5 , UPdGa 5 , and UPtGa 5 . On the other hand, for 
M=Co, Rh, Ir, Fe, Ru, and Os, both Xa and Xc are al- 
most independent of temperature, since these are Pauli 
paramagnets, but in these cases Xc is somewhat larger 
than Xa- Nonetheless, we again observe a tendency for 
the anisotropy to become larger for the progression from 
3d to 4d and then to hd M ions. Thus, it is characteris- 
tic of U-115 compounds that the magnetic properties are 
sensitive to the choice of transition metal ions. 

Furthermore, UGa3, which is the mother compound of 
UMGas, has also provided intriguing experimental re- 
sults. It has been reported that UGa3 exhibits a G- 
type AF metallic phase in the low-temperature region^i 
but a "hidden" ordering different from the magnetic 



one has been suggested by resonant X-ray scattering 
measurements^ Unfortunately, orbital ordering in UGa3 
is not yet confirmed experimentally, but it may be an 
interesting possibility to understand the result of reso- 
nant X-ray scattering experiment on UGa3 based on the 
orbital-ordering scenario. 

Here we note that one must pay close attention to the 
meanings of "spin" and "orbital" in /-electron systems. 
Since they are tightly coupled with each other through 
a strong spin-orbit interaction, distinguishing them is 
not straightforward in comparison with d-electron sys- 
tems. This point can create serious problems when we 
attempt to understand microscopic aspects of magnetism 
and superconductivity in /-electron compounds. Thus, 
it is necessary to carefully define the terms "orbital" and 
"spin" for / electrons in a microscopic discussion of mag- 
netism and superconductivity in uranium compounds. 
In order to overcome this difficulty, we have proposed 
to employ a j-j coupling scheme to discuss /-electron 
systems 4^ 

Let us note the advantages of the j-j coupling scheme 
at the outset. First, it is quite convenient for the in- 
clusion of many-body effects using standard quantum- 
field theoretical techniques, since individual /-electron 
states are clearly defined. In contrast, in the LS cou- 
pling scheme we cannot use such standard techniques, 
since Wick's theorem does not hold. Second we can, in 
principle, include the effects of valence fluctuations. In 
some uranium compounds, the valence of the uranium 
ion is neither definitely U 3+ nor U 4+ , indicating that the 
/-electron number takes a value between 2 and 3. In 
the j-j coupling scheme this is simply regarded as the 
average number of / electron per uranium ion. 

In this paper, then, we propose a new scenario based 
on the j-j coupling scheme in order to understand the 
magnetic properties of uranium compounds, of which 
UGa3 and UMGas are considered to be typical exam- 
ples. A microscopic model constructed from the j-j cou- 



pling scheme is analyzed using an unbiased method, such 
as the Lanczos technique. For UGa3 we obtain a phase 
diagram, including the paramagnetic (PM) state and sev- 
eral types of AF phases. In particular, we note that the 
PM state is adjacent to an AF metallic phase, consistent 
with the experimental results for uranium compounds 
having the AuCu3-type crystal structure. We also dis- 
cuss the orbital structure of the AF metallic phase. By 
taking account of the effects of two-dimensionality, we 
find that the change in magnetic structure from G- to 
A-type AF phases is reproduced, consistent with exper- 
imental results for UMGas. Moreover, it is shown that 
the observed trends of the magnetic anisotropy within 
the UMGa5 series can be understood from the present 
results. 

The organization of this paper is as follows. In Sec. II, 
we will introduce a microscopic model Hamiltonian for- 
mulated on the basis of the j-j coupling scheme. In 
Sec. Ill, numerical results for UGa3 and UMGas are dis- 
cussed in detail and compared with experimental results. 
In Sec. IV, future developments are discussed and the 
paper is summarized. In the Appendix, we analyze the 
local /-electron configuration in order to examine the 
validity of the j-j coupling picture vis-a-vis experimen- 
tal results on the level scheme in certain /-electron ma- 
terials. Throughout the paper, we use units such that 
H=kB = l- 



II. MODEL HAMILTONIAN 

In this section we derive a microscopic Hamiltonian for 
UGa3 and UMGas based on the j-j coupling scheme. Al- 
though it is difficult to determine the valence of the ura- 
nium ion, here we assume that the formal valence is U + , 
including three / electrons per ion. By considering the 
crystalline electric field (CEF) potential and Coulomb in- 
teractions, we then assign three electrons to states in the 
j=5/2 sextet. 
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FIG. 1: Configurations of / electrons in the j-j coupling 
scheme for (a) Celn3 and (b) UGa3. Note that up and down 
arrows denote pseudospin states in order to distinguish the 
two states of the Kramers doublet. 

and excited states, respectively, with an energy difference 
of 12meVii£ Thus, we take T7 to be lower for the present 
considerations, as shown in Fig. 1(a). 

In the j-j coupling scheme for UGa3, we accommodate 
three electrons in the one-electron energy states and 
IV There are two possibilities, i.e., low- and high-spin 
states, depending on the Hund's rule interaction and the 
splitting between the and Tg levels. Noting that the 
effective Hund's rule interaction is small in the j-j cou- 
pling scheme^ the low-spin state should be realized, as 
shown in Fig. 1(b). In fact, this low-spin state is also 
consistent with the LS coupling scheme. Details regard- 
ing electron configurations in the j-j coupling scheme are 
discussed in the Appendix. 

In the electron configuration shown in Fig. 1(b), the 
T7 level is fully occupied to form a singlet. If this 
level is located well below the r 8 , the occupying electrons 
will not contribute to the magnetic properties. Thus, we 
may be allowed to ignore the T7 electrons for our present 
purposes, but this simplification will be examined more 
carefully in the following subsection. 



B. Suppression of T7 



A. Local /-electron states 

First, we define the one /-electron states in the AuCu3- 
type cubic crystal structure. The effects of tetragonality 
will be discussed later. From the work of Hutchings for 
the case of cubic symmetry^ we identify two eigen en- 
ergies as -2405^ for the T 7 doublet and 120B® for the 
Tg quartet, where B® is a cubic CEF parameter. 

In order to proceed with the discussion, it is necessary 
to know which is lower, T7 or Tg, in the one /-electron 
picture. For some crystal structures it is possible to de- 
termine the level scheme from intuitive discussions of /- 
electron wavefunctions and the positions of ligand ions. 
However, this is not the case for the AuCu3-type crystal 
structure. For this case we invoke certain experimental 
results on Celn3, a typical AuCu3-type Ce-based com- 
pound, where T-j and Tg, have been reported as ground 



First, we discuss the /-electron kinetic term given by 
#kin = - t* T ,ff Ta f i+aTa , (1) 

i,R,T,r' ,cr 

where f- lTIJ is the annihilation operator for an /-electron 
with pseudospin a in the r-orbital at site i, and t* T , is 
the /-electron hopping matrix element between r-and r'- 
orbitals along the a direction. Indices a and b denote the 
r| and Tg orbitals, respectively, while c indicates the T^. 
Here, the pseudospin has been introduced to distinguish 
the two states of the Kramers doublet. 

When we evaluate t^. T , for nearest-neighbor hopping 
via the a bond, it is given by 

/ 3/4 -n/3/4 0\ 
% T , = t -V3/4 1/4 , (2) 
V / 
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for the cc-direction, 

/ 3/4 V3/4 0\ 
t\ T , = t V3/4 1/4 , (3) 

\ o oo/ 

for the y-dircction, and 

/ \ 

t*r> =0^0, (4) 
\ / 

for the z-direction. Note that t z =t in the cubic system. 
In this paper, t is taken as the unit of energy. Later, 
t z will be treated as an extra parameter in order to take 
account of the tetragonality. We immediately notice that 
non-zero hoppings occur only among Tg orbitals. Since 
the 1?7 orbital has nodes along the cubic axes, it is lo- 
calized in the present tight-binding approximation. This 
further justifies the suppression of T?. 

One may question the foregoing approach on the 
grounds that it is based entirely on nearest-neighbor hop- 
ping processes. In order to address this concern, it is 
useful to examine the results of band-structure calcu- 
lations, e.g., for Celn 3 (Ref. HI) and UGa 3 (Ref. 
Note that both results have been obtained assuming the 
system is in the paramagnetic state. In order to focus 
on the / electron components of the energy band, we 
concentrate on the bands around the T point near the 
Fermi level. For Celn3, the energy band dominated by 
character is found to be lower than the Fg-dominated 
band, consistent with the local level scheme in Fig. 1(a). 
An important point is that the Fermi level intersects the 
r7-dominant band, indicating that the Fermi surface is 
mainly composed of Ti electrons hybridized with Ga-ion 
p electrons. 

On the other hand, for UGa3, the band is also lower 
than the Tg band, but here the Fermi level crosses the Tg 
band. Thus, the T7 band appears to be fully occupied, 
consistent with the j-j coupling level scheme, as shown 
in Fig. 1(b). Since the main contribution to the Fermi 
surface comes from Tg electrons, it is natural to dwell on 
the Tg bands and ignore the occupied T7 bands in giving 
further consideration to many-body effects. 

It is useful to consider the Fermi-surface structure of 
-ffkin in comparison with that of the band calculations. In 
Fig. 2, we show Fermi-surface sheets derived from -ffkin 
for the case of (n)=l, where (n) indicates the average 
number of Tg electrons per site. Owing to the multi- 
orbital nature of the problem, we observe two Fermi- 
surface sheets. A cube-like Fermi surface is centered on 
the r point and a second sheet is composed of three con- 
nected tubes along the three cubic axes. 

In the band-structure calculation results, two Fermi 
surface sheets are also observed^ One sheet surrounds 
the r point with a small hole at the center. Another is a 
large sphere-like Fermi surface centered at the R point. 
Since the carrier number is simply fixed at unity per site 
in the tight-binding model, it is difficult to obtain perfect 




FIG. 2: Fermi-surface sheets of the Tg tight-binding model 
for {n) — l. Note that we obtain two Fermi-surface sheets, (a) 
and (b), due to the multi-orbital nature of the model. Note 
also that both sheets are depicted in the electron picture. The 
bounding box indicates the first Brillouin zone for a simple 
cubic lattice. The F point is located at the center of the box, 
while the apices denote R points. 

agreement. Except for details, however, -ffkin can repro- 
duce the Fermi-surface structure of the band calculation 
in spite of its simplification. 

C. Ts model 

Following the above discussion, the effective Hamil- 
tonian for uranium compounds will be expressed by a 
rs-orbital degenerate Hubbard model, given as 

H = -ffkin + -ffcEF + -ffint, (5) 

where -ffkin denotes the Tg part of Eq. (1) and -ffcEF is 
written as 

H C ef = -A y^(nj g - nib). (6) 

i 

Here n iT = ^ CT n iTcr and n iT(T = ff Ta .fi T<r . The level split- 
ting A between Tg orbitals is introduced in order to 
consider certain tetragonal CEF effects for the case of 
UMGa5. When we analyze the magnetic properties of 
UGa3, this term is not needed. For the time being, A is 
then set to be zero. 

The Coulomb interaction term H ln t is given by 

-ffint = U ^2 ™iTT«-irl + V ^ n ia n ib 
i.r i 

+ J/% f \ T a f{ T ' a' fir a' fir' a 
i,(T,er' r^r' 

i t^t' 

where U, U' , J, and J' denote intra-orbital, inter-orbital, 
exchange, and pair-hopping interactions, respectively. 
These are expressed in terms of Racah parameters, and 
the relation U=U'+J+J' holds due to the rotational in- 
variance in orbital spaced For (i-electron systems, one 
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also has the relation J=J'. When the electronic wave- 
function is real, this relation is easily demonstrated from 
the definition of the Coulomb integral. However, in the 
j-j coupling scheme the wavefunction is complex, and J 
is not equal to J' in general. For simplicity, we shall 
assume here that J= J', noting that essential results are 
not affected. Since double occupancy of the same orbital 
is suppressed owing to the large value of U, pair-hopping 
processes are irrelevant in the present case. 

III. RESULTS 

Among several possible methods to analyze the present 
microscopic model, in this paper we have employed the 
exact diagonalization technique. Although this has the 
drawback that it is difficult to enlarge the system size, 
this method offers the clear advantage that it is possible 
to deduce the magnetic structure in an unbiased manner. 
In order to discuss the ground-state properties, it is useful 
to calculate both the spin and orbital correlations, which 
are defined, respectively, by 

S{ q ) = {l/N)Y J {°l°!V q - {i - i \ (8) 

i j 

with <t?= J2 T ( n ^ - n- 1T {)/2 and 

T{q) = {l/N)Y J {r-;r?)e^\ (9) 

with T? = {rii a — nib)/2. Here N is the number of sites. 

A. Cubic system 

First let us consider the cubic system. Due to se- 
vere limitations of computer memory, our calculations 
are restricted to a 2x2x2 cube. Nonetheless, essen- 
tial points on the spin structure can be captured. In 
Fig. 3(a), as a typical result for spin correlation, we show 
S(q) as a function of J for U'=6. From the changes 
in the dominant component, we can define three regions 
as (I) J<0.45, (II) 0.45<J<2.05, and (III) J>2.05. In 
region I, the dominant component of S(q) appears at 
q=(ir, 7r, 7r), indicating a G-type AF structure. Defini- 
tions of spin structure are shown in Fig. 3(b). Note 
that for C-type [q=(n, ir, 0), (ir, n,0), (it, 7r,0)] and A- 
type [q=(n, 0,0), (0,7r,0), (0,0, 7r)] structure, S(q) has 
the same value for each q due to the cubic symmetry. 
In region II, 0.45< J<2.05, the cubic symmetry seems to 
be broken. Of course, this is spurious because of the 
smallness of the system, but for reasons given below, this 
phase is considered to be "metallic" , and it is convention- 
ally called the PM phase. Note that the AF correlation 
with q — (tt,it,tt) is still dominant for 0.45<J<0.8, in- 
dicating that the system is in a metallic magnetic phase. 
For 0.8< J<2.05, A-type AF correlations turn out to be 
dominant. 
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FIG. 3: (a) Spin correlation S(q) as a function of J for U'—6. 
(b) Spin structures in the FM, A-type AF, Gtype AF, and 
G-type AF phases. 



Here we discuss the evidence for a crossover between 
insulating and metallic behavior given by the spin corre- 
lation function. In an insulating phase, the spin struc- 
ture is essentially determined by a round trip for an elec- 
tron just between neighboring sites, leading to orbital- 
dependent superexchange interactions. Thus, cubic sym- 
metry is maintained in the spin correlation function even 
for a small-sized cluster such as a 2x2x2 cube. However, 
in the PM metallic phase, electrons tend to gain kinetic 
energy by moving around the whole system. Such a mo- 
tion depends sensitively on the anisotropy of the shape 
of the orbital. Specifically, in a small-sized cluster, the 
spin correlation function in the PM phase depends on the 
choice of the basis set for orbitals, and the cubic symme- 
try of the spin correlations appear to be broken owing to 
the smallness of the system. Thus, the spurious violation 
of cubic symmetry in the spin correlations is considered 
to be a signal for a metallic phase. Note that the metal- 
insulator boundary itself depends on the system size. In 
the thermodynamic limit, cubic symmetry should exist in 
the metallic phase, since the effect of orbital anisotropy is 
smeared by averaging over the whole system. In fact, the 
band dispersion relation has cubic symmetry, as is easily 
checked by diagonalizing in momentum space. 

When we further increase the value of J, the spin cor- 
relations recover cubic symmetry in region III, J>2.05. 
Again we find an insulating phase, but the spin structure 
is now considered to be A-type AF. The change from G- 
to A-type AF phase with increasing J can be understood 
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FIG. 4: (a) Orbital correlation as a function of J for U = 6 
Meanings of the symbols are the same as those in Fig. 3(a) 
Orbital arrangements (b) for J<0.7 and (c) for J>0.7. 



in terms of the competition between kinetic and magnetic 
energies, in analogy with manganites^ Namely, for large 
J there is an energy gain for an FM spin pair on neigh- 
boring sites, but not for an AF spin pair. Thus, there is 
a tendency for the occurrence of ferromagnetism in the 
large-J region. When we compare the magnetic energy 
of nearest-neighbor pairs, we notice that the number of 
FM pairs in the A-type AF phase is larger than for the 
G-type phase. Thus, for large values of J the A-type 
phase appears. 

In Fig. 4(a) we show typical results for orbital cor- 
relations as a function of J for U'—6. Correspond- 
ing to the changes in the dominant component of the 
spin correlation, we again see three regions. In region 
I (J<0.45), T(q) has two degenerate dominant compo- 
nents of <7=(7r, 0,0) and (0, 7t, 0). The orbital pattern is 
given by a mixture of two A-type orbital correlations, as 
schematically shown in Fig. 4(b). Since the Hamiltonian 
is not invariant under orbital exchange, T(q) should not 
have the same magnitude for q=(Tr, 0,0), (0, 7r,0), and 
(0,0, 7r). On the other hand, in region III (J>2.05) the 
(•7T, 7r, 7r) component is dominant, suggesting a G-type or- 
bital pattern as shown in Fig. 4(c). 

Now we focus our attention to the region II, in which 
we can see the crossover in the orbital pattern between 
Figs. 4(b) and (c). Note that the crossover point is dif- 
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FIG. 5: (a) Spin correlation as a function of U' for J=0. 
(b) Orbital correlation as a function of U' for J—0. (c) Phase 
diagram for UGa3 obtained by the exact diagonalization. The 
region of J>U' is ignored, since it is unphysical. See Fig. 3(b) 
for the definitions of abbreviations. Here "PM-G" indicates 
the PM phase with enhanced (ir, n, tv) spin correlation. 



ferent from that between G- and A-type AF phases ob- 
served in S(q). This may be related to the appearance 
of two distinct transition temperatures of UGa3 around 
at Xn-— Namely, one is related to ordering of spin and 
another is originating from orbital degree of freedom. An 
interesting point is that both orbital patterns, Figs. 4(b) 
and (c), are antiferro-like in the basal plane. This may 
be detected in future by careful experiments. 
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In order to consider the spin and orbital structure of 
uranium compounds, let us focus on the small- J region. 
For J=0, we show S(q) as a function of U' in Fig. 5(a). 
In the PM region for small U' such as for U'<2, there is 
no dominant component in S(q), but for £7'>2, (tt,tt,tt) 
correlation gradually grows. The PM phase with dom- 
inant (-7T, 7T, 7r) spin correlation is defined as "PM-G" in 
this paper. With further increasing U' , eventually the 
ground state becomes G-type AF insulating for U'>A. 

In the experimental results for UX3, when the lattice 
constant becomes large in the order of X=Si, Ge, Sn, and 
Pb, the system is changed from Pauli paramagnetic to AF 
metallic^ For X=A1, Ga, In, and Tl, a similar change 
has been reportedi22i2i If we simply consider that the 
increase of the lattice constant leads to the decrease in the 
effective hopping, U' /t becomes large with the increase 
of the lattice constant. Then, the present phase diagram 
in the small- J region is consistent with the experimental 
results for UX3. 

In Fig. 5(b), the orbital correlation as a function of 
U' is shown For J=0. In the region corresponding to 
PM-G, both orbital patterns, Figs. 4(b) and (c), are pos- 
sible in principle, but we note that the pattern (b) ap- 
pears for larger value of V in the PM-G region near 
the G-AF insulating phase. Since UGa3 is an AF metal, 
it may be natural to consider that the orbital correla- 
tion corresponding to Fig. 4(b) can be detected in actual 
compounds. Note also that two types of orbital arrange- 
ment, Figs. 4(b) and (c), should be distinguished by the 
structure along the c-axis. This point may be clarified in 
future experiments. 

After we have performed calculations for several pa- 
rameter sets, the ground-state phase diagram is com- 
pleted on the (U',J) plane, as shown in Fig. 5(c). The 
PM phase exists for large parameter space and in the 
boundary region between PM and G-type AF states, we 
can see the PM phase with dominant (ir,Tr,ir) spin cor- 
relation. Note that such a PM-G region is not specific to 
the case of J=0, since it appears even when we increase 
the Hund's rule interaction. 

Here we briefly discuss the phases in the large- J region. 
We observe an interesting similarity with the phase dia- 
gram for undoped manganites RMn03^ in which mobile 
e g -electrons are tightly coupled with the Jahn- Teller dis- 
tortions and the background t 2g spins. Note that the 
present Hamiltonian is just equal to the e g electron part 
of the model for manganitesii^ In the so-called double- 
exchange system with large Hund's rule coupling between 
e g and t 2g electrons, the Jahn- Teller distortion suppresses 
the probability of double occupancy and it plays a simi- 
lar role as the interorbital Coulomb interaction U'. The 
AF coupling among t 2g spins, Jaf, controls the FM ten- 
dency in the e g -electron phases. Roughly speaking, large 
(small) Jaf denotes small (large) J. Then, we see an 
interesting similarity between Fig. 5(c) and the phase di- 
agram for manganites, except for the PM region. Espe- 
cially, a chain of the transition, FM— >A-AF— >C-AF— >G- 
AF, occurs with decreasing J (increasing Jaf)- Again 



we stress that the present Tg model for /-electron sys- 
tems is essentially the same as the e g orbital model in 
the d-electron systems. It is interesting to observe com- 
mon phenomena concerning orbital degree of freedom in 
/-electron systems. 



B. Tetragonal system 

In the previous subsection, we have analyzed our model 
for the case of a cubic system in an effort to understand 
UGa3. Our analysis has yielded a magnetic, metallic 
phase with antiferro-like orbital correlations, consistent 
with the experimental results. In order to extend this dis- 
cussion to tetragonal systems such as UMGa5, it is nec- 
essary to introduce two new ingredients into the model 
Hamiltonian. 

One is a non-zero value for A, which is a level splitting 
between two of the orbitals. Under a tetragonal CEF, 
there are two levels and one Tg. Of these the Tq state 
is the same as Tg for the cubic system. The two states 
are given by linear combinations of the J z =±3/2 and 
=F5/2 states, which can also be expressed as an admixture 
of T7 and Tg. For simplicity, we introduce A as a splitting 
energy between the Tg, orbitals. 

Secondly, it is necessary to change the hopping am- 
plitude along the z-axis. In UMGas, the MGa2 layer is 
sandwiched between two UGa3 sheets, suggesting that 
the hopping amplitude of /-electrons along the z-axis 
should be reduced from that in UGa3. It is difficult to 
estimate this reduction quantitatively, since one must in- 
clude correctly the hybridization with d-electrons in the 
transition metal ions and with p-electrons in the Ga ions. 
Thus, we manage the anticipated reduction by simply 
treating t z as a parameter. 

In Figs. 6(a) and (b), typical results for spin and or- 
bital correlations are shown for ^=0.8, J=0, and U'=3.5. 
For |A|^>1, a G-AF phase is observed, since the Hamil- 
tonian is effectively reduced to a single-band model at 
half-filling, and the superexchange interaction stabilizes 
the G-AF state. However, the A-AF phase appears 
for — 0.39<A<— 0.17 near the orbitally degenerate re- 
gion. The mechanism of the appearance of the A-AF 
phase in the negative A region will be discussed later, 
in terms of an orbital-based scenario. Regarding the or- 
bital structure, for |A|^>1, simple fcrro-orbital (FO) or- 
dered phases are obtained. In the narrow region where 
— 0.23<A<— 0.05, we have identified an antiferro orbital 
(AFO) pattern, as shown in Fig. 4(c). Although the 
crossover point from an FO to an AFO pattern deviates 
slightly from the A-AF and G-AF phase boundary, it is 
basically considered that the A-AF phase appears in the 
region with an FO pattern. 

In Fig. 6(c), we show the phase diagram in the (A,t z ) 
plane for J=0 and U'—3.5, in which the ground state for 
A=0 and t z =l is magnetic metallic, as seen in Fig. 5(c). 
It is found that an A-type AF phase appears in the neg- 
ative A region for i 2 >0.68. Note that the appearance of 
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FIG. 6: Spin (a) and orbital (b) correlation functions vs. A 
for t z —Q.S, with J—0 and f/'=3.5. (c) Phase diagram of the 
magnetic structure in the (A, t z ) plane for J—0 and V = 3.5. 
(d) Ferro orbital pattern in the A-type AF phase. 



the A-AF phase is not sensitive to t z as long as ?i 2 >0.68. 
Rather, A seems to play a key role in controlling the 
change of the magnetic phase. Here we recall the exper- 
imental fact that UNiGa5 exhibits a G-type AF phase, 
while UPtGas shows an A-type^ Thus, it is necessary to 
relate the effect of A to the difference in magnetic struc- 
ture found between UNiGas and UPtGa 5 . Although t z 
may differ among U-115 compounds, we focus here on 
the effect of A. 

Let us now discuss the reasons for the appearance of 
an A-AF phase. For negative values of A, we easily ob- 



tain FO ordering composed of Tg orbitals, as illustrated 
in Fig. 6(d). For electrons to gain kinetic energy of mo- 
tion along the z-axis, it is necessary to place the AF spin 
arrangement along this same axis. In the FM spin con- 
figuration, electrons cannot move along the z-axis due to 
the Pauli principle, since hopping occurs only between 
r| orbitals along the z-axis. On the other hand, in the 
xy plane o-orbital electrons can hop to neighboring a- 
orbitals with a significant amplitude, which is larger than 
that between neighboring 6-orbitals. Thus, in order to 
gain kinetic energy, electrons tend to occupy a-orbitals 
even in the FO state composed of &-orbitals, as long as | A| 
is not so large. When we explicitly include the effects of 
the Hund's rule interaction J, electron spins should have 
FM alignment between neighboring sites in order to gain 
energy in hopping processes from b- to a-orbitals. Con- 
sequently, a FM spin configuration is favored in the xy 
plane. In fact, in spite of the FO state for — 0.39<A<0, 
we can see a significant component of T(n, it, tt). In cases 
with antiferro orbital correlations, spin correlation tends 
in general to be FM, as has been widely recognized in 
orbitally degenerate systems. 

Here we mention a relation of A to the magnetic 
anisotropy in U-115 materials. For UPtGas with the A- 
AF phase, Xa is larger than Xc whereas this anisotropy 
is not pronounced in UNiGas with the G-AF phase. 9 
An analysis for the high-temperature region based on 
LS coupling yields the J z =±l/2 Kramers doublet as 
the ground state among the dectet of J— 9/2 (L=6 and 
5=3/2) ^ xhe states with J z =±l/2 in the LS coupling 
scheme have significant overlap with // C J0) and 

fiblfLlfbl\°) in the 3-3 coupling scheme. Accordingly, 
by the present definition A should be negative to place T\ 
below rjj. If the absolute value of A(<0) becomes large, 
F| is well separated from and the magnetic anisotropy 
will consequently become large. Thus, a change from G- 
to A-type AF phase is consistent with the trends of mag- 
netic anisotropy in UNiGas and UPtGas. 

Finally, we make a brief comment about the effect of t z . 
Following the above discussion, the A-AF phase should 
appear even for small t z . However, in the persent cal- 
culation it disappears for i z <0.68, a critical value which 
seems to be rather large. Such a quantitative point de- 
pends on the system size, and we note that it is necessary 
to perform the calculation in the thermodynamic limit. 
This is a problem for future consideration. 



IV. DISCUSSION AND SUMMARY 

In this paper, we have proposed a calculational model 
with active orbital degrees of freedom in an effort to un- 
derstand the magnetic structure of uranium compounds 
from a microscopic point of view. In order to construct 
such a model, we have incorporated the j-j coupling 
scheme, in which one-electron states are defined first, and 
Coulomb interactions are included subsequently. This 
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approach is consistent with the itinerant picture for / 
electrons. By using an exact diagonalization technique, 
we have found a magnetic metallic phase with antiferro- 
like orbital correlations for UGa3 and a change in the 
magnetic structure from G- to A-type AF phases for 
UMGa 5 . 

In an effort to understand the magnetism and super- 
conductivity of /-electron systems from a microscopic 
standpoint, we have carried out a study of an orbitally 
degenerate model. While such investigations are just be- 
ginning, we already see a number of opportunities for 
future work along this path. Concerning issues directly 
related to the present paper, it is highly recommended 
that calculations be carried out in the thermodynamic 
limit, in order to confirm the present exact diagonaliza- 
tion results. For instance, the magnetic susceptibility 
should be evaluated in the random phase approximation 
or fluctuation-exchange method. With such an approach, 
the magnetic structure can be discussed by detecting the 
divergence in the magnetic susceptibility. This is one of 
our future tasks. Another problem is how to establish 
the effective reduction of t z in considering the case of 
UMGa5. In such systems, MGa2 sheets are interspersed 
between UGa3 layers, but the main process may occur 
through the Ga ions>i£ To analyse this, it is necessary to 
treat a three-dimensional f-p model with explicit consid- 
eration of U and Ga ions. This is another problem for 
future investigation. 

We comment briefly on the CEF levels of 115 com- 
pounds. In this paper, we have considered the UMGa5 
systems by introducing a splitting energy A between 
Tg orbitals. This procedure did not include changes in 
the basis wavefunctions. Nonetheless, in following the 
present strategy, a better way to consider U-115 might 
be to accommodate three electrons in the CEF levels of 
Ce-115. Regarding the latter systems, there are some in- 
teresting results in the literature. In analyses of the sus- 
ceptibility and specific heat^iS 5 - for example, the ground 

(2) 

level has been considered to be I7 , consistent with 
Xa<Xc m Ce-115, and the first excited level has been sug- 
gested to be Ty 1 ^, not r 6 . Regarding CeRhIn 5 , neutron 
scattering experiments have been performed in spite of 
the inclusion of indium^ identifying energy levels in the 
sequence lf>, lf\ and r 6 , in agreement with previous 
analyses^ If we simply put three electrons into the CEF 
levels of CeRhlns, assuming Hund's rule coupling to be 
small, it may be difficult to understand the appearance 
of A-type AF phases, since A becomes positive within 
the level scheme of Ce-115. 

A simple way to understand such a discrepancy is to 
remark that Ce-115 contains In and M = Co, Ir, and Rh, 
while the U-115 considered here contains Ga and M = 
Ni, Pd, and Pt. Different ions lead to different effects 
on the level scheme. In fact, the magnetic anisotropy 
of UCoGa5 (Pauli paramagnet) is small compared with 
that of UNiGas, and \ c is slightly larger than Xa- In 
UFeGa5, which is also a Pauli paramagnet, we have found 
that Xa<Xc, opposite to what is observed in UNiGa5 



and UPtGas. We envisage that in UMGas, A is posi- 
tive for M=Fe, slightly positive for M=Co, and negative 
for M=Ni. Moreover, |A| increases in the sequence 3d, 
Ad, and 5d transition metal ions, consistent with the mag- 
netic anisotropy. 

In summary, we have analyzed an orbitally degenerate 
model appropriate for UGa3 and UMGas using an ex- 
act diagonalization technique. What we have found is a 
magnetic metallic phase with antiferro-likc orbital corre- 
lations. By introducing tetragonal effects such as level 
splitting and a reduced hopping amplitude along z-axis, 
we have reproduced the change in the spin structure from 
G- to A-type AF phases, corresponding to UNiGa5 and 
UPtGa 5 . 
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APPENDIX A: LEVEL SCHEME 

In this Appendix, we examine a local /-electron state 
based on the j-j coupling scheme, comparing it with that 
obtained from the LS coupling scheme and from exper- 
imental results. Concerning the level scheme for the f 1 
electron state, we consider two cases. 

First, we consider the AuCu3-type cubic crystal struc- 
ture with one / electron per site. As already men- 
tioned in the main text, a typical material is Celn3, in 
which T7 and Tg are the ground and first excited states, 
respectively, 15 as shown in Fig. 7(a). If we accommodate 
one more electron to consider the / 2 configuration, imme- 
diately there appear two possibilities, "low" and "high"- 
spin states. When the CEF splitting energy between Ti 
and Tg levels is smaller than the Hund's rule coupling, 
the second electron should be accommodated in the Tg 
levels. In the situation in which one is in the T>? and the 
other in the T%, & T4 or T§ triplet appears for the / 2 
state in general, but under special conditions, a T3 dou- 
blet can occur. On the other hand, if the CEF splitting 
is larger than the Hund's rule interaction, then the / 2 
ground state is formed from two Ty electrons, leading to 
a Ti state. When we compare this Ti state with that 
in the LS coupling scheme, we notice that it is given by 
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FIG. 7: Electron configurations in the j-j coupling scheme for 
AuCu3-type compounds, (a) Celn3, (b) Prln3, and (c) Ndln3. 



FIG. 8: Electron configurations in the j-j coupling scheme for 
rare-earth hexaborides, (a) CeB6, (b) PrBs, and (c) NdB6. 



a mixture of J=0 and J— A states, but the J— A compo- 
nent is found to be dominant. Note also that Ti is the 
antisymmetric representation of T7 x FV. 

Since we do not know the exact value of the Hund's 
rule interaction in /-electron compounds, it is difficult to 
determine the f 2 state by purely theoretical arguments. 
In this case, we have to refer to data on actual materi- 
als. Fortunately, we have the example of Prln3, a typical 
f 2 material with AuCu3-type crystal structure. From 
several experimental results, Ti has been confirmed to 
be the ground level in Prln3^i Thus, the low-spin state 
should be taken for the AuCu3-type structure in the j-j 
coupling scheme. 

Here the reader may pose a naive question: Is the 
Hund's rule interaction really that small in /-electron 
systems? We have already had to answer this question in 
a previous paper ^ In a word, we are considering the ef- 
fective Hund's rule interaction in the j-j coupling scheme, 
not in the LS coupling scheme. The original form for the 
Hund's rule interaction is written as — Jh« 2 , where s\ de- 
notes the operator for the "real" /-electron spin at site 
i and Jjj is the Hund's rule interaction among / elec- 
trons in £=3 orbitals. In order to transfer this to the j-j 
coupling scheme, it is convenient to use the well-known 
relation Sj=(gj— l)ji, where gj is the Lande's <?-factor 
and ji is the operator for the total angular momentum 
of j=5/2 at site i. From standard textbooks, we easily 
obtain gj=6/7, indicating that Si=— (l/7)ji. Thus, the 
Hund's rule term in the j-j coupling scheme is rewritten 
as — Jcffjf with J c ff=Jn/49. Note then that the magni- 
tude of the Hund's rule interaction is effectively reduced 
by the factor 1/49 in the j-j coupling scheme. Even 
if Jn=leV, Jeff is reduced to be about 200K, which is 
comparable with the CEF splitting energy. Thus, it is 
possible to have the low-spin state in the j-j couping 
scheme. 

Next, we take a further step to the / 3 state by adding 

one more / electron. Since T-j is fully occupied to form 

Ti, the next electron should be placed in the r 8 state as 

shown in Fig. 7(c), clearly indicating that there exists an 

active orbital degree of freedom. The / 3 state composed 

(2) 

of two and one Tg electron is expressed as T y s ' in the 
terminology of group theory. When we again consider ac- 
tual materials, Ndln3 is found to be a typical / 3 material 



with the AuCu3-type crystal structure. In experiments, 
it has been established that Ig is the ground level^ as 
we have found with the present j-j coupling scheme. 

Let us turn our attention to another crystal structure 
in which Tg is lower than T^ in the f 1 configuration. Typ- 
ical materials are the rare-earth hexaborides RBg with 
R=Ce, Pr, and Nd. As is well known, the ground level 
of CeBg is rg, indicating that the quadrupolar degree of 
freedom plays an active role in this material^ In fact, 
anomalous behavior related to quadrupolar ordering has 
been suggested by several exeprimental results. 

First, we note that the level splitting between Tg and 
is assumed to be larger than the Hund's rule interac- 
tion. When we accommodate two electrons in Tg orbitals, 
the triplet (r 5 ), doublet (r 3 ), and singlet (ri) states are 
allowed. Among these, owing to the effect of the Hund's 
rule interaction, even if it is small, the T$ triplet should 
be the ground state. This has actually been observed 
in PrB6^2*^i Further, in order to consider NdB6, yet an- 
other electron is put into the Tg orbital, making a total of 
three. Alternatively, we may say that there is one hole in 

the Tg orbital. Such a state is found, again, to be char- 
(2) 

actcrized by T 8 . Experimental results on NdB@ have 
actually been reported which lead to the ground state of 
p( 2 ) , 30,32 Thus, when r 8 is the ground state for the one 

/-electron case, we obtain T§ for the f 2 and F g 2 ^ for the 
/ 3 configurations. 

We have shown that the ground states deduced from 
the j-j couping scheme are consistent with experimental 
results. However, in order to explain the experimental re- 
sults quantitatively, it is unavoidable to analyze the CEF 
levels using the LS coupling scheme. What we would like 
to stress is that even in a localized system, the symme- 
try of the ground level can be understood via the j-j 
coupling scheme. We need to recognize the limitations of 
the j-j coupling scheme when we treat a local electronic 
state. For instance, to consider the / 3 state, we simply 
put three electrons into the CEF level scheme which is 
determined with the f 1 configuration. Thus, the wave- 
function of the / 3 state is uniquely determined. How- 
ever, in an actual situation, the dectet labelled by J=9/2 
{L—Q and 5=3/2) is split into two Tg and one Tq orbital. 
The ground-state wavefunctions will then depend on the 
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two CEF parameters and B§£& As mentioned above, 
in order to explain experimental results on localized /- 
electron materials, ones should analyze the CEF effects 
in the system using the LS coupling scheme. In this pa- 
per, however, the electronic states are considered with 
an itinerant picture based on the j-j coupling scheme. 
Thus, it is important to check that the local electronic 
state formed by / electrons in this way is consistent with 
the symmetry of the state obtained with the LS coupling 
scheme. 



In summary, it has been shown that the ground states 
of the f 2 and / 3 configurations can be qualitatively re- 
produced by accommodating / electrons in the CEF lev- 
els of a corresponding f 1 material, provided that the CEF 
level splitting is larger than the Hund's rule interaction. 
Thus, the j-j coupling scheme works even in the localized 
case. Accordingly, we believe that a microscopic theory 
can be developed in which we discuss the magnetism and 
superconductivity of /-electron compounds in terms of 
the j-j coupling scheme. 
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